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APPARATUS AND METHOD FOR CORRECTING ERRORS GENERATED BY 
A LASER WITH NON-IDEAL TUNING CHARACTERISTICS 

CROSS-REFERENCE TO RELATED APPLICATIONS 
This application claims the benefit of U.S. Provisional Patent Application Serial 
No. 60/459,856, titled, "Apparatus and Method for Correction Errors Generated by a 
Laser with Non-Ideal Tuning Characteristics," filed April 2, 2003, and is hereby 
incorporated by reference in its entirety. 

FIELD OF THE INVENTION 
The present invention relates to lasers having non-ideal tuning characteristics and 
their use in optical measurements; primarily measurements using Optical Frequency 
Domain Reflectometry. In particular, it relates to an apparatus and-method for correcting 
errors generated by lasers having non-ideal tuning characteristics. 

BACKGROUND OF THE INVENTION 
U.S. 5,798,521 to Froggatt, discloses an apparatus and method for measuring 
strain of gratings written into an optical fiber. By measuring the complete spectral 
response of the Bragg grating, the strain at each point in the grating can be measured. 
The apparatus comprises a control system, a data acquisition (DAQ) circuit, detectors, a 
laser controller, and a coherent light source or tunable laser. The reference and 
measurement fringes are detected and sampled such that each sampled value of the 
reference and measurement fringes is associated with a corresponding sample number. 
The wavelength change of the laser, for each sample number, is determined by 
processing the signal from the reference optical fiber interferometer. Each determined 
wavelength change is matched with a corresponding sampled value of each measurement 
fringe. Each sampled measurement fringe of each wavelength sweep is transformed into 
a spatial domain waveform. The spatial domain waveforms are summed to form a 
summation spatial domain waveform that is used to determine location of each grating 
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with respect to a reference reflector. A portion of each spatial domain waveform that 
corresponds to a particular grating is determined and transformed into a corresponding 
frequency spectrum representation. The strain on the grating at each wavelength of 
optical radiation is determined by determining the difference between the current 
wavelength and an earlier, zero-strain wavelength measurement. For this application, the 
measurement of the instantaneous wavelength of the tunable laser by the single reference 
interferometer is sufficiently accurate. 

U.S. patent 6,376,830 to Froggatt et al. is directed toward a system and method 
for measuring the transfer function of a guided wave device. In particular, the NxN 
scalar transfer function elements for an N-port guided wave device are measured. Optical 
energy of a selected wavelength is generated at a source and directed along N reference 
optical paths having N reference path lengths. Each reference optical path terminates in 
one of N detectors such that N reference signals are produced at the N detectors. The 
reference signals are indicative of amplitude, phase and frequency of the optical energy 
carried along the N reference optical paths. The optical energy from the source is also 
directed to the N-ports of the guided wave device and then on to each of the N detectors 
such that N measurement optical paths are defined between the source and each of the N 
detectors. A portion of the optical energy is modified in terms of at least one of the 
amplitude and phase to produce N modified signals at each of the N detectors. At each of 
the N detectors, each of the N modified signals is combined with a corresponding one of 
the N reference signals to produce corresponding N combined signals at each of the N 
detectors. A total of N 2 measurement signals are generated by the N detectors. Each of 
the N measurement signals is sampled at a wave number increment Ak so that N 
sampled signals are produced. The NxN transfer function elements are generated using 
the N 2 sampled signals. Reference and measurement path length constraints are defined 
such that the N combined signals at each of the N detectors are spatially separated from 
one another in the time domain. 

Because U.S. 6.376,830 is directed toward optical instrumentation, the accuracy 
requirements on the phase measurements are quite severe. The measurement is based on 
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the assumption that a long interferometer can be used to precisely monitor the 
instantaneous wavelength of a tunable laser. This assumption is valid under most 
regimes of operation, and to the degree-of-accuracy generally required. However, when 
lasers with significant tuning-speed variations are used coupled with long (>20m) paths 
and a requirement of phase accuracies on the order of milli-radians, then, the assumption 
is no longer valid (as is the case with optical instrumentation). The breakdown of the 
assumption occurs when the timescale of the tuning speed variation occurs on the time 
scale of the delay in the interferometer. In turn, it becomes necessary to correct for errors 
generated by the optical source. 

U.S. 6,566,648 to Froggatt describes an apparatus and method for measuring 
strain of gratings written into an optical fiber. Optical radiation is transmitted over one or 
more contiguous predetermined wavelength ranges into a reference optical fiber network 
and an optical fiber network under test to produce a plurality of reference interference 
fringes and measurement interference fringes, respectively. The reference and 
measurement fringes are detected, and the reference fringes trigger the sampling of the 
measurement fringes. This results in the measurement fringes being sampled at 2n 
increments of the reference fringes. Each sampled measurement fringe of each 
wavelength sweep is transformed into a spatial domain waveform. The spatial domain 
waveforms are summed to form a summation spatial domain waveform that is used to 
determine location of each grating with respect to a reference reflector. A portion of each 
spatial domain waveform that corresponds to a particular grating is determined and 
transformed into a corresponding frequency spectrum representation. The strain on the 
grating at each wavelength of optical radiation is determined by determining the 
difference between the current wavelength and an earlier, zero-strain wavelength 
measurement. The apparatus and method disclosed herein fails to disclose an auxiliary 
interferometer that corrects for residual errors resulting from the laser. 

An object of the present invention is to provide an apparatus and method for 
correcting for errors generated by a laser with non-ideal tuning characteristics. 
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Another object of the present invention is to provide an apparatus and method for 
correcting errors generated by a laser with non-ideal tuning characteristics that employs 
at least one auxiliary interferometer in parallel with a sampling interferometer and a 
measurement interferometer. 

5 

SUMMARY OF THE INVENTION 
By the present invention, an apparatus for correcting errors generated by a laser 
with non-ideal tuning characteristics is presented. For the purpose of the instant 
application and the appended claims, non-ideal is defined as a laser having a tuning speed 
10 that is not sufficiently linear to provide phase measurements that meet the user's 

requirements. 

The apparatus comprises a laser having non-ideal tuning characteristics. At least 
three interferometers are positioned in an operable relationship to the laser. A sampling 
interferometer and at least one auxiliary interferometer correct for residual errors 

15 resulting from the laser. A measurement interferometer makes a measurement. A signal 

acquisition system is positioned in an operable relationship to each interferometer. A 
processor is positioned in an operable relationship to the signal acquisition system. 

A method for correcting errors generated by a laser with non-ideal tuning 
characteristics is also presented. In practicing the method, a laser having non-ideal 

20 tuning characteristics is provided. At least three interferometers are positioned in an 

operable relationship to the laser. A sampling interferometer and at least one auxiliary 
interferometer correct for residual errors resulting from the laser. A measurement 
interferometer makes a measurement. A signal acquisition system is positioned in an 
operable relationship to each interferometer and a processor is positioned in an operable 

25 relationship to the signal acquisition system. 

Additional objects and advantages of the invention will be set forth in part in the 
description which follows, and in part will be obvious from the description, or may be 
learned by practice of the invention. The objects and advantages of the invention will be 
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obtained by means of instrumentalities in combinations particularly pointed out in the 
appended claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The accompanying drawings illustrate a complete embodiment of the invention 
according to the best modes so far devised for the practical application of the principles 
thereof, and in which: 

FIG. 1 is a schematic representation of an embodiment of the invention in its 
simplest form where three interferometers are employed. 

FIG. 2 is a schematic representation of the invention where multiple 
interferometers are used to measure difference functions as continuous variables of time. 

FIG. 3 is a schematic representation of the invention where one of the correcting 
interferometers-is constructed in a Michelson interferometer configuration using Faraday - 
Rotating Mirrors (FRM's). 

FIG. 4 is a schematic representation of the invention in a Michelson 
interferometer configuration as it is incorporated into a system. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Fiber optic interferometers are often used to measure various changes within a 
sample such as pressure, temperature, strain and stress. When conducting a 
measurement, light is launched through the optical fiber by a laser. The light passes 
through an optical Device Under Test (DUT) and an amplitude and phase versus time 
signal is generated by the interferometer. If the tuning speed of the laser is not 
sufficiently linear to provide phase measurements that meet the user's requirements, the 
laser is deemed to have non-ideal tuning characteristics. 

The present invention lends itself to the correction of errors generated by the 
tuning variations that exist in a tunable laser. Referring now to the figures, where similar 
elements are numbered the same throughout, in general, the invention is an apparatus for 
correcting errors generated by a laser with non-ideal tuning characteristics. As shown in 
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FIG. 1, the apparatus comprises a laser 10 having non-ideal tuning characteristics. Any 
laser known to those of ordinary skill in the art may be used in the present invention. 
Preferably, the laser is tunable and mode-hop free. 

At least three interferometers 20, 30, 40 are positioned in an operable relationship 
to the laser. A sampling interferometer 20 and at least one auxiliary interferometer 30 
correct for residual errors resulting from the laser. In practice, the sampling 
interferometer 20 corrects for most of the errors caused during the laser sweep. The 
auxiliary interferometer 30 corrects for errors not captured by the sampling 
interferometer 20. A measurement interferometer 40 is used to measure a sample. The 
figure depicts a preferred embodiment where the measurement interferometer 40 has a 
device under test positioned therein. 

Each interferometer 20, 30, 40 generates an analog signal versus time. Signals 
from 30 and 40 are collected through a signal acquisition system 50 with each sample - 
triggered by the signal from 20, such as a data acquisition device, which is positioned in 
an operable relationship to each interferometer. A National Instruments PCI-61 15 12-bit, 
10 Megasamples per second data acquisition card was used in this embodiment. A 
processor 60 is positioned in an operable relationship to the signal acquisition system 
50. The processor 60 is also positioned in an operable relationship to the laser 10 to 
signal the laser when to begin the test. 

FIG. 2 depicts an alternative arrangement of the invention wherein two auxiliary 
interferometers 30, 35 are employed in combination with the sampling interferometer and 
the measurement interferometer. When it is desirable to increase precision, additional 
auxiliary interferometers are added. As each additional auxiliary interferometer is added, 
additional residual errors resulting from the laser are corrected. Thus, the more 
corrections, the higher the precision of the measurement to be made. 

In U.S. Patent No. 6,376,830 a system and method for measuring the transfer 
function associated with single-port guided wave devices or the transfer function matrix 
of a multi-port guided wave device, e.g., Bragg gratings, couplers, etc. was proposed. 
The problem with this system and method is that it is based on the assumption that a 
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single interferometer is adequate to correct for tuning speed variations in the laser. 
Moreover, such a system employed interferometers to make measurements from a device 
having light going through it. However, this system fails to sufficiently eliminate errors 
that are generated by the source of the light. Therefore, to solve this problem, begin with 
5 a presumed laser field of: 

i/f(t) = e i ( a) ° t+at2+acos<a \ 

where, a cos cot represents the tuning speed error as time harmonic phase error. 
10 The laser frequency increases linearly with a small oscillation of amplitude, a , and 

frequency, co . The interference fringes formed by delaying a portion of this field and 
adding it to an undelayed version is then given by: 

p(t) = y/(ty (0+ r(t - x y (f - r)+ 2\^(t]\^(t - r) cos(^ c ) 

15 

where r is the length of the optical delay in time. The phase of the fringe term is then 
given by: 



"fringe 

20 



fringe = + ^ + Q C0S 0)1 ~ ^0 (* ~ T )~ a (* ~" T Y ™ ° C0S 0)(t ~ t) . 

Collecting terms. 



& fringe = 0) Q T - (XT 2 + 2<XTt + tf(cOS Q)t - COS Co{t - t)) 

Trigonometrical ly expanding the cosine of a difference. 

25 

$ fringe ~ a) o T ~ ar2 + ^ art + # ( cos ~ cos ®t cos cot - sin cot sin cot) 
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Collecting terms again 



$fiin& ^^r-ar 1 + 2aU + a[cos#tf(l- cos cor)- sin cot sin cor] 



By dropping the constant phase term (non time-varying) in front, the following equation 
is obtained: 



This fringe signal is used to sample another fringe signal from an interferometer with 
some length that is unequal to the sampling interferometer. The presumption follows that 
this is the phase of the sampled measurement interferometer and will be directly 
proportional to the phase of the sampling interferometer. Let the sampling interferometer 
have a delay, r v 3 and measurement interferometer have a phase, r m . The two phases are 
then given by: 

fa = 2ar J + a[cos cot(\ - cos cor s ) - sin cot sin cor s ] 



The corrected phase error, as caused by sampling fa as a function of fa is then given by, 



<P fringe = 2(XTt + 4 C0S ^K 1 ~ 



cos<5>r 



)- sin cot sin a>r] 



and 



(f> m = 2arJ + a [cos cot(l - cos cor m )- sin cot sin cor m ] . 



A =6 - T -^d> 

rem rm r.v 



where a linear relationship is assumed such that it is proportional to the delays. 



Substituting in: 
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And collecting terms. 



= a 



cos cot{\ - cos (or m ) - — cos cot(l - cos 6)T S ) - sin cot sin co r m + — sin cot sin cor s 

T T 

K V 



The frequency of the laser tuning speed variation, co , multiplied by the interferometer 
delays, r v and r OT is assumed to be a small but non-zero number. The sine and cosine 
terms are expanded to second order such that 
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= a 



cos cot 



2 2 \ / 

# r I r f 
1-1 + =- — ^cos*tf 



v 



1-1 + 



V 



- ^yr m sin cot + —cor s sin 



15 



Collecting terms 



f 2 2 2 2 \ 

2 r. 2 



COS Ct)t 



And finally 



A-=^-[(«"i-V,)cosfljf] 



20 



This expression results in a phase amplitude error which takes the form of a parabola 
with a maximum value at / 2 r x and zeros at 0 and r, . This maximum value is: 



9 



Docket No.: LT 182-1 



Patent Application 



r max g 

Putting this in more standard units, 

/ W f 124.7 ^2^(lxlO- 9 GHzmzVr)r; 124.7 A , x 2 

- w t" = l^r^J — (ads = WtfJ&f* h 

where, r , is the average tuning rate of the laser in nm/s, Ar is the amplitude of the 
deviation of the laser tuning rate in nm/s, and A/l is the period of the variation. Using the 
Agilent 81640A laser, the tuning rate is 40 nm/s, the period of variation is 0.03 nm 

AtfL = 7 l lV \ M40X400) 2 = 3.3rad 

^ max (8xl0 9 )(0.030) V A ; 

Correction for the induced error can be made if another interferometer is constructed with 
a known single-path difference. From this interferometer, a measurement of the phase 
error due to tuning-speed variations can be obtained. 

^=^[fe-^.v)cos^] 

The ratio of the phase errors is then given by: 

6 T 2 —T T 

rem m m s 

6 T 2 —T T 

Tea a * a \ 

Therefore, if a known interferometer is constructed and the phase error incurred 
when sampling this delay is measured, the phase error for all delays can be predicted. 
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Generally, electronic delays in the detection circuit and signal conditioning can add 
delays and render the above equation an approximation. The ratio of the auxiliary phase 
to the measurement phase error can be found experimentally for any given apparatus by 
measuring the phase errors directly for a variety of measurement lengths. When this is 
done, it has been found that the variation of the ratio with length is quadratic as in the 
above equation. 

The second order error due to laser tuning speed fluctuations can then be removed 
by using the following formula: 



t — r t 

"m, corrected rm 2 ' a 

r — T T 

a as 



where $ m and % , are the phases of the signals measured using the sampling - ~ 
interferometer to trigger the sampling of these signals. 

Considering the expression for the magnitude of the phase error due to tuning 
variation: 



one may observe that it depends quadratically on the frequency of the tuning error. As a 
result, slow errors (such as those in the low audio regions — 100 of Hz) do not contribute 
significantly. These errors can be removed by unwrapping the auxiliary phase signal, 
curve fitting it, and subtracting off the fitted curve. 

Referring now to FIG. 2, the labeled path lengths (ao, aj, bo, bi, co,ci, and do ,di) 
are assumed to incorporate the full path from the source to the detector (not shown). If 
the signals w, x, y, and z are assumed to be phase signals, then the interferometers 
measure difference functions as continuous variables of time, such that: 
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w(t) = f(t-a ( )-f(t-a 0 ) 
x(t)=f{t-b x )-f{t-b 0 ) 
yit)=f{t-c x )-f{t-c 0 ) 
z (t) = f(t- d] )-f(t-d 0 ) 

It is assumed that w(t) is the trigger interferometer phase, and z(t) is the measurement 
interferometer phase. Therefore, with no device present it is expected that the phase of 
z(t) can be placed at a constant value with the proper corrections. 

The initial estimate of the proper correction q 0 (t) is given by: 

Provided that: 
f(t) = r 0 +r l t + r 2 t 2 , 

then q 0 (t) is a constant: 

q 0 (t) = C 

and independent of t . This amounts to the constant tuning rate approximation. 
Alternatively, by allowing a tuning rate that changes, another order can be added: 



/(')=r 0 +r/ + r/+r 3 f\ 

and, if all of the reference paths are equal, a 0 
q(t) that is first order in / . 

^ 0 (0 = C + 3r 3 (^-fl 0 Xrf I -a ! >. 



-b Q =c 0 =d 0 , then an error is picked up on 



12 



Docket No.: LT182-1 



Patent Application 



This error can be eliminated by recognizing that another term dependent upon t is 
obtained by evaluating: 

a } -a 0 

a linear term in / results: 
p o (0 = C + 3r 3 (6 1 -a o X6 I -a I > 

This signal is used to correct (to first order) for tuning speed variations in the laser, 
providing a new phase estimator: 

a x -a 0 

Given the balanced reference paths described above, this correction eliminates any errors 
due to third order phase effects in the laser (i.e. changes in the scan rate). If, however, 
there are fourth order terms, then another linear term is added. That linear term has a 
coefficient: 

q x if) = C + 6r 3 lVf°^ * ^ + b * d » + d * a ° ~ a « b > ~ a ^~ b » d > * 

{b l -b 0 )-{a l -a 0 ) 

If the fourth order term: 
f{t) = r 0 +r t t + r/ + r/ + r/ 

is present and the reference paths are balanced, then the error term is, 



a \~ a Q 



(d x -d 0 \d,-a x ) 
(b,-b Q \b,-a x ) 
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and the error goes to zero when the measurement length, d } , is equal to the reference 
length (i.e. zero length difference), when the measurement length is equal to the trigger 
delay, a x , and when the measurement length is equal to the correcting (auxiliary) 
interferometer length. 

The remaining interferometer signal, y(t) , is used to further correct for this fourth 
order term: 

q, (r ) = z(t) - fa~4Xfr-4Xq-4) ^ _ fa-4Xq-4Xq-4) ^ _ fa-4Xq-<0fa-4) ^ 

fa -cfifa -bja, -cj fa -hfa -<hfa fa> -<hfa -<0fe "4) 

This process is repeated, adding more interferometers to obtain arbitrarily good 
estimations. 

FIG. 3 depicts another embodiment of the invention where the sampling 
interferometer 20 is constructed in a Michelson interferometer configuration using 
Faraday Rotator Mirrors (FRM's) 70. Faraday Rotator Mirrors 70 ensure that full 
interference is achieved over a wide range of tuning and all input polarization states. The 
first order interferometer is generally the longest, and is thus most susceptible to effects 
from optical polarization. The auxiliary interferometers 30, 35 coupled with the sampling 
interferometer 20 correct for residual errors resulting from the laser tuning speed 
variations. The measurement interferometer 40 is shown in a preferred embodiment 
where it has a device under test positioned therein. 

FIG. 4 depicts an apparatus for correcting errors generated by a laser with non- 
ideal tuning characteristics where the sampling interferometer 20 is a Michelson 
interferometer configuration using Faraday Rotating Mirrors (FRM's) 20. An auxiliary 
interferometer 30 and the sampling interferometer 20 correct for residual errors resulting 
from the laser 10. The measurement interferometer 40 is shown in a preferred 
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embodiment where it has a device under test (DUT) positioned therein. A signal 
acquisition system 50 is positioned in an operable relationship to each interferometer 20, 
30 5 40 such that it receives the analog phase versus time signal generated by each 
interferometer. A processor 60 is positioned in an operable relationship to the signal 
acquisition system. The processor 60 uses phase information to correct for residual errors 
generated by the laser. The first order correction is carried out by triggering data 
acquisition points as described in U.S. Patent No. 6,376,830 to Froggatt et al. and in U.S. 
Patent No. 6,566,648 to Froggatt, the specifications of each are hereby incorporated by 
reference in their entirety. The calculations of the phase signals from the acquired analog 
signals are carried out by the methods descried in U.S. Patent No. 5,798,521 to Froggatt 
and is hereby incorporated by reference in its entirety. 

In practicing the method of the present invention, a laser having non-ideal tuning 
- characteristics is provided. At least three interferometers are positioned in an operable- 
relationship to the laser. A sampling interferometer and at least one auxiliary 
interferometer correct for residual errors resulting from the laser, where a measurement 
interferometer makes a measurement. In a preferred embodiment, the measurement 
interferometer has a device under test positioned therein. A signal acquisition system is 
positioned in an operable relationship to each interferometer. A processor is positioned 
in an operable relationship to the signal acquisition system and uses phase information 
from the sampling and auxiliary interferometers to correct for measurement errors 
generated by the measurement interferometer. In an alternative embodiment, when it is 
desirable to further correct for errors, at least one more interferometer is employed. 

When errors occur due to third order phase effects in the laser, such errors are 
corrected by the equation: 

a,-a 0 

where w, x and z are phase signals and where: 



d, -d n 



{d x -d 0 \d x -a x ) 
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*{/)=. f{t-a x )-f(t-a a ) 
x(t) = f(t-b t )-f(t-b 0 ) 
z{t) = f(t-d i )-f(t-d 0 ) 

and ai, ao , bi, bo ,and di, d 0 are path lengths. 
Fourth order terms are corrected by the equation: 

q (t) - z(t ) ^° - d[ ^ - d ' - dl ! u(f ) ^° - d> ^' - d > ^ - d ^ x{t) ^° ~ d< ^' - d ' ^ ~ ^ ^ >(f ) 

2 (^-0|X«i-^X«|-c,) (a 0 -*,X*.- a iXA- c i) (oo-^X c r-^X c .-^) 

where w, x, y and z are phase signals and where: 
w{t) = f{t- a] )-f{t-a 0 ) 

*W=/M)-/('-*b) 

X/) = /(f-c,)-/(f-c 0 ) 
z(/)=/(/-J,)-/(/-J 0 ) 

and ai, ao ? bi ? b 0 ,cj, co,and di, d 0 are path lengths. 

The above description and drawings are only illustrative of preferred 
embodiments which achieve the objects, features and advantages of the present invention, 
and it is not intended that the present invention be limited thereto. Any modification of 
the present invention which comes within the spirit and scope of the following claims is 
considered part of the present invention. 



16 



